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Abstract 

A perturbative approach for arbitrary choices of the equation of state of the universe is introduced 
in order to treat scenarios for mass varying neutrinos (MaVaN's) coupled to the dark sector. The 
generahzed criterion for the apphcabihty of such an approach is expressed through a constraint on 
the coefficient of the linear perturbation on the dark sector scalar field. This coefficient depends on 
the ratio between the variation of the neutrino energy and the scalar field potential. Upon certain 
conditions, the usual stationary condition found in the context of MaVaN models together with 
the perturbative contribution can be employed to predict the dynamical evolution of the neutrino 
mass. Our results clearly indicate that the positiveness of the squared speed of sound of the coupled 
fluid and the model stability are not conditioned by the stationary condition. 
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The ultimate nature of the dark energy and its connection with the Standard Model (SM) 
states is one of the most intriguing issues related with the negative pressure component 
required to understand the accelerated expansion of the universe. If from one hand, the 
most natural candidate to couple a SM singlet quintessence-like scalar field is the Higgs 
sector , an exciting alternative is the coupling to neutrinos as in the context of the 
so-called mass varying neutrino (MaVaN) models |3|-l7|. This last possibility is particularly 
interesting since the coupling of neutrinos to the dark energy scalar field component may lead 
to a number of significant phenomenological consequences. Actually, this coupling renders 
the MaVaN mechanisms fairly natural. Indeed, if the neutrino mass mj, is generated by 
the dynamical value of a cosmologically active scalar field instead of through a vacuum 
expectation value (VEV) it would be an evolving quantity. 

For the purpose of analyzing the minimal energy density contribution of neutrinos to the 
cosmic budget, a stationary condition (SC) has been devised so to establish a relationship 
between and the scalar field potential. In this work we consider the possibility of neu- 
trino masses arising from an interaction with the real scalar field that drives the accelerated 
expansion of the universe. In order to obtain a model along these lines we consider the Ma- 
VaN contribution to the energy conservation equation of the cosmic fluid as a perturbation 
that iteratively modifies the background fluid equation of state. This approach is taylored 
to satisfy the phenomenological requirement that the contribution of neutrinos to the energy 
density evolution is subdominant in the matter dominated era. The scalar field is modified 
by a linear perturbation due to the neutrino mass coupling, which is turned on when the 
neutrinos become non-relativistic (NR). In such a proposal (see Ref. |3] for a detailed discus- 
sion), one is compelled to abandon the SC even though this condition can be reconciliated 
with our perturbative approach through a particular neutrino mass dependence on the scale 
factor. 

Another relevant related issue is that dark matter is most often not considered in the 
formulation of the MaVaN models. However, the possibility of treating dark energy and 
dark matter in a unified scheme naturally offers this possibility. The generalized Chaplygin 
gas (GCG) is particularly relevant in this respect as it is shown to be consistent with 
known observational constraints, from CMB [10], supernova 
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12], gravitational lensing 



surveys [13], gamma ray bursts [14] and cosmic topology [l5]. Our analysis considers the 
coupling of MaVaN's to the underlying scalar field of the generalized Chaplygin gas (GCG) 
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model. Since the neutrino contribution is perturbative, we obtain a small deviation from 
the stability condition characteristic of the unperturbed (GCG) equation of state. 

Given that cosmological neutrinos have not been observed so far, hints at about their 
contribution to the energy density are relevant, albeit somewhat theoretical. The neutrino 
energy density and pressure are expressed through a Fermi-Dirac distribution function with- 
out a chemical potential, /(<?), where q = T^o being the neutrino background temperature 
at present. Thus, 



-1/2 



where we have assumed a flat FRW cosmology, that the scale factor of the universe at present 
is = 1 (the sub-index denotes present-day values) and, for simplicity, we considered only 
one non-vanishing neutrino mass. By observing that 

my{4,)^^^^^-^ = {pu{a,(t>) - 3]9^(a, (/))), (2) 

and from Eq. ([21), one can obtain the energy-momentum conservation for the neutrino fluid 

Puia, 0) + 3H{p^{a, (/)) + pu{a, 0)) = 0— 7- ^ — , (3) 

where H = a/a is the expansion rate of the universe and the overdot denotes differentiation 
with respect to cosmic time (' = d/dt). 

It is important to emphasize that the coupling between cosmological neutrinos and the 
scalar field as specified in E g. (|2l) is restricted to times when neutrinos are NR, i. e. 
^dm^{(j)) ~ ^^^""^ ^ — !■ other hand, as long as neutrinos are relativis- 

tic {Ty{a) = Tya/a » my{(t){a))) , the decoupled fluids should evolve adiabatically since the 
strength of the coupling is suppressed by the relativistic increase of pressure {p^ ~ "ipu)- In 
this case, one should have 

p^ + 3i7(p<^ + p^) =0, (4) 
and a similar equation for neutrinos. For the scalar field background fluid one has 

P0 = y + V{4>), 

= Y - V[c^). (5) 



Treating the system of NR neutrinos and the scalar field as a unified fiuid (UF) which 
adiabatically expands with energy density pup = Pu + P<t> and pressure pup = Pu + P<f) leads 
to 

puF + ^H{puF + Puf) = =^ P0 + 3if (p0 + = (6) 

where the last step is derived from the substitution of Eq. ([3]) in the first equation in ([6]). 

It is well known that the relative contribution of the energy densities components of the 
universe with respect to the one of the dark energy sector is on its own a p roblem. The 
assumptions proposed in Ref. j4| and subsequently developed elsewhere [gI. IigIiS] introduce 
a SC which allows circumventing the coincidence problem for cosmological neutrinos, by 
imposing that the dark energy is always diluted at the same rate as the neutrino fiuid, that 
is, 

dV{(j)) _ dm^ dpu , s 

dcj) d(f) drrii, 

This condition introduces a constraint on the neutrino mass since it promotes it into a 
dynamical quantity, as indicated in Eq. ([6]). In this context, the main feature of the FNW 
scenario is that, in what concerns to dark sector, it is equivalent to a cosmological 
constant-like equation of state and an energy density that is as a function of the neutrino 

n 

mass [7|]. As already mentioned, the effectiveness of this coupling is limited by values of the 
scale factor greater than ajv^, where a^vB parameterizes the transition between the relativistic 
and NR regimes. 

The assumption of a universe with the dark energy governed by the equation of state 
Pff) = —pti) implies, through Eq. ([7]), that Ph = V , and allows to recover the SC, exactly 
as obtained in Ref. Q]. It is clear that the relevance and the considerations about the 
constraints on the equation of state are conditioned by the underlying scalar field potentials 
and their adequacy concerning neutrino mass generation. In fact, once one assumes that 
p^ = —Pfj,, the neutrino mass evolution and the form of the potential become directly 
entangled by the SC. We believe that an alternative approach for treating deviations from 
the SC is needed since the constraint that it ensues, namely P4, + Pcf, = 0, is quite restrictive 
as it disregards the kinetic energy contribution of the scalar field. Although this contribution 
is at present negligible in comparison to the one of the potential, it can dominate and play 
a role in the evolution of the universe at earlier times. 

Given these considerations, instead of beginning from Eq. ([6]), we use the unperturbed 
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Eq. dl]), and establish the conditions for treating the neutrino couphng in a perturbative 
way. To exemphfy this suggestion more concretely, let us consider the unperturbed equation 
of motion for the scalar field 

+ 3if0+^ = O, (8) 

and assume that the effect of the coupling of the neutrino fluid to the scalar field fluid is 
quantified by a linear perturbation ecf) (|e| << 1) such that — )• ^ (H-e)0. It then follows 
a novel equation for energy conservation 

(^ + 3i70+ = (9) 
dip dip dniy 

After some straightforward manipulations one can easily find that 

^"^ + (^'^)^' 
which, upon substitution into Eq. ([9]), and use of Eq. (|8]), leads to 



dW{(j)) dV{(l)) 



dm^ dpu drriy . 



dp drriu 

where the perturbative character of the neutrino mass term is assumed when setting the 
last approximation in the above equation. Finally, we obtain for the value of the coefficient 
of the perturbation 

drriv dpi, 



d(f> dm I, 



d fldV{<j)) 



(12) 



which for consistency is required to satisfy the condition |e| << 1. Upon fulfilling all known 
phenomenological requirements, the above result allows us to address a wide class of scalar 
field potentials and related equations of state for various candidates for the dark sector 
(dark energy and dark matter), which through the SC would be incompatible with realistic 
neutrino mass generation models. 

Notice that the equation for the conservation of energy, Eq. ([9]), can be reobtained by 
simply redefining p^p as 

Pu. = 1<P' + Kff (13) 
so that the usual definition of an effective potential V-es in terms of = ^^^^ + 

dmy dpu 
dij) dm.1, 



inbL^Pji. ig now valid for any value of 
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In order to verify under which conditions Eq. ([7]) agrees with our perturbative approach 
for a given background equation of state, the coefficient of the hnear perturbation should be 
given by 

(|e|«l). (14) 



,2A (^i^lM 

d<j> \(f> dip 

This means that we must search for a neutrino mass dependence on the scale factor for 
which the above condition is satisfied. Thus, once one sets the equation of state for the dark 
sector, there will be a period at late times for which the SC and the perturbative approach 
match. In particular, this feature can be reproduced by the GCG equation of state given 
by, 

p=-^.po(^y, (15) 

where < < 1 and < a < 1 js, 9|. Inserting the above equation into the unperturbed 
energy conservation Eq. (jlj), one can integrate it to obtain 



P<f> = Po 



'1 - A., 



i/(i+a) 



(16) 



Q3(i+a) 

from which P(^(a) is easily computed, and by Eqs. ([5]), the scalar potential V{4'{a)) [ll| is 
written as 



A, + 



:i - A.. 



-a/{l + a) 



2 A 



A, 



a 



3(l+a) 



(17) 



a3(i+a) 

Clearly, given a potential, the explicit dependence of m(0) on a can be immediately obtained 
from Eq. ([7]). Furthermore, it is necessary to determine for which values of the scale factor 
the neutrino-scalar field coupling becomes important. For convenience we set the value of 
a = affR for which Pi,^nr = Pu,ur holds, usually established by the condition of > T^, that 
parameterizes the transition between the NR and the ultra-relativistic (UR) regime. In fact, 
this takes place when 

T 

mu{a) = X—, (18) 

a 

where x is a numerical factor estimated to be about % ~ 94 considering that 
Pu/Pcvit = mo [eV^]/(94 /i^ [ey]), where h is the value of the Hubble constant in terms of 
100 km Mpc~^ . Such a correspondence between a^vA and mo is illustrated in the Fig. [T] 
for the particular value of a = 1/2. Considering the whole set of parameters that character- 
ize the background fiuid, one notices that is rather difficult to see that the maximal value 
assumed by the e parameter corresponds to its present-day value, as seen in the Fig. [2l 



Mass -mo 
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FIG. 1: Present-day values of the neutrino mass ttiq and the corresponding values of Ojvii for which 
the transition between the NR and UR regimes takes place in the GCG phenomenological scenario 
with Q = 1/2 and variable Ag. The increasing graylevel corresponds to increasing values of nio, for 
which we have marked the boundary values for mo = 0.05 eV, 0.1 eV, 0.5 eV, 1 eV, 5 eV. 

For a more accurate comparison between our perturbative approach and the SC frame- 
work, it would be interesting to express analytically the bound (e < 1) for the coupling 
function in terms of the model parameters. Considering dark sector scalar fields and poten- 
tials with a simple analytical dependence on the scale factor, we can compare the constraint 
imposed on the potential by the breaking of the perturbative approach with the known 
limits derived in the literature for stable MaVaN scenarios jo]. However, due to its some- 
what complex analytical description, the GCG fluid is not the most suitable for this kind of 
quantitative comparison. 

We observe that the interval of parameters A^, rrio and eventually a, for which our 



approximation can be applied (e < 1), is valid for a > a^^n and severally constrainec 



by the 



imposition af^^ < 1. For phenomenologically viable values of (0.7 < < 1) [10|, ll2| one 
finds that 0.1 < e < 1. In general terms, just under quite special circumstances the usual SC 
and the perturbative contribution of MaVaN's match. In the original MaVaN scenario j4|, 
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FIG. 2: Maximal value of the linear perturbation coefficient e as a function of the scale factor and 
of the GCG parameters = 0.7 and a = 1/2. For illustration purposes we have set mo = 0.5 eV. 
The increasing gray level corresponds to increasing values of |e|, for which we have marked the 
boundary values for |e[ = 0.01, 0.1, 1. 

the SC corresponds to the adiabatic solution (if^ ^ d^V/dcp"^) of the scalar field equation 
of motion. In this case, the kinetic energy terms of the scalar field can be safely neglected. 
The consistency of our perturbative scenario with the stationary condition can be achieved 
only when the kinetic energy contribution is not relevant at late times. From Eqs. f lT6|) . f lT7|) 
one can easily recover the kinetic energy component of the GCG fluid in terms of ^(^(a)) 
and p^, 



which can be set equal to (1 — A^)/2 at present times (a = 1). Turning back to Eq. flTBl) . 
one notices that, for A, = 0, the GCG behaves always as matter, whereas for A, = 1, it 
behaves always as a cosmological constant. Consequently, it is natural that the relevance of 
the kinetic energy term at present times is suppressed when the parameter A^ gets close to 
unity, which enlarges the agreement between our perturbative approach and the SC analysis. 
In what concerns the results of our approach, we illustrate the increasing the neutrino 
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FIG. 3: Model independent perturbative modification on squared speed of sound as a function 
of the scale factor for the neutrino-GCG coupled fluid in comparison with the adiabatic GCG fluid 
for = 0.7 with a = 1, 1/2, 1/3, 10"", for a present-day value of neutrino mass, mo = 0.5 eF. 

mass with the scale factor for a set of phenomenologically consistent parameters for the 
GCG in Fig. |3l Interestingly, from Fig. [3l for mo = 0.5 eV^, a fairly typical value, we can 
see that stable MaVaN perturbations correspond to a well defined effective squared speed 
of sound, 



of sound, up to the limit where the perturbative approach breaks down. However, we find 
that as far as the perturbative approach is concerned, our model does not run into stability 
problems in the NR neutrino regime. In opposition, in the SC treatment, where neutrinos 
are just coupled to dark energy, cosmic expansion in combination with gravitational drag 
due to cold dark matter have a major impact on the stability of MaVaN models. Usually, for 
a general fluid for which the equation of state is known, the dominant behaviour on c\ arises 
from the dark sector component and not by the neutrino component. For the models where 
the SC (cf. Eq. ([7])) implies in a cosmological constant-type equation of state, p,^ = — p,/,, 
one obtains = — 1 from the very beginning of the analysis. 




> 0. 



(20) 



The greater are the mo values, the more important are the corrections to the squared speed 



9 



Moreover, our perturbative approach is in agreement with the assumption that the cou- 
phng between neutrinos and dark energy (and/or dark matter) is weak. We have observed 
that the stabihty condition related to the squared speed of sound of the coupled fluid is 
predominantly governed by dark energy equation of state. Such a troublesome behaviour 
should in fact have already been observed in the context of the FNW scenario as the SC 
implies that cl = —1 from the very start and the role of recovering the c^f^^^^^^ > condition 
is relegated to the neutrino contribution 18|]. The loosening of the stationary constraint 
Eq. ([7]) emerges from the dynamical dependence on (p, more concretely due to a kinetic 
energy component The knowledge of the background fluid equation of state for the dark 
sector (the GCG in our example), and the criterion for the applicability of the perturbative 
approach, do allow us to overcome the negative problem, independently from the neutrino 
mass dependence set by the SC. 

As a general final remark, we mention that coupling a dark sector field (dark energy or 
dark energy plus dark matter) to SM states may bring important insights about the physics 
beyond the SM. Neutrino cosmology, in particular, is a fascinating example from where 
salient questions concerning SM particle phenomenology can be resolved or, at least, better 
understood if one considers the coupling to the dark sector scalar field. 



Acknowledgments 

A. E. B. thanks the financial support from the FAPESP (Brazilian Agency) grant 
07/53108-2 and the hospitality of the Physics Department of the Instituto Superior Tecnico, 
Lisboa, Portugal. O. B. would like to acknowledge the partial support of Fundagao para 
Ciencia e Tecnologia (Portuguese Agency) under the project POCI/FIS/56093/2004. 



[1] B. Patt and F. Wilczek, l arXiv:hep-ph/0605188t 

[2] O. Bertolami and R. Rosenfeld. lar Xiv:n7n8.l"784l [hep-ph]. 

[3] P. Gu, X. Wang and X. Zhang, Phys. Rev. D68, 087301 (2003). 

[4] R. Fardon, A. E. Nelson and N. Weiner, JCAP 0410 005 (2004). 

[5] J. Lesgourgues and S. Pastor, Phys. Rept. 429, 307 (2006). 

[6] O. E. Bjaelde et al, JCAP 0801, 026 (2008). 

10 



[7] A. E. Bernardini and O. Bertolami. larXiv:0712. 15341 [astro-ph]. 

[8] A. Kamenshchik, U. Moschella and V. Pasquier, Phys. Lett. B511, 265 (2001); N. Bilic, G. 

B. Tupper and R. D. Viollier, Phys. Lett. B535, 17 (2002). 
[9] M. C. Bento, O. Bertolami and A. A. Sen, Phys. Rev. D66, 043507 (2002); 
[10] M. C. Bento, O. Bertolami and A. A. Sen, Phys. Rev. D67, 063003 (2003); Phys. Lett. B575, 

172 (2003); L. Amendola, F. Finelh, C. Burigana and D. Carturan, JCAP 0307, 005 (2003). 
[11] O. Bertolami, A. A. Sen, S. Sen and P. T. Silva, Mon. Not. Roy. Astron. Soc. 353, 329 (2004). 
[12] M. C. Bento, O. Bertolami, A. A. Sen and N. C. Santos, Phys. Rev.D71, 063501 (2005). 
[13] P. T. Silva and O. Bertolami, Ap. J. 599, 829 (2003); A. Dev., D. Jain and J. S. Alcaniz, 

Astron. Astrophys. 417, 847 (2004). 
[14] O. Bertolami and P. T. Silva, Mon. Not. Roy. Astron. Soc. 365, 1149 (2006). 
[15] M. C. Bento, O. Bertolami, M. J. Rebougas and P. T. Silva, Phys. Rev. D73, 043504 (2006). 
[16] R. D. Peccei, Phys. Rev. D71, 023527 (2005). 

[17] A. W. Brookfield, C. van de Bruck, D. F. Mota and D. Tocchini-Valentini, Phys. Rev. D73, 
083515 (2006). 

[18] R. Takahashi and M. Tanimoto, JHEP 0605, 021 (2006). 



11 



